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Abstract. Our former study on the empirical relations between the Fourier parameters of the light curves of the 
fundamental mode RR Lyrae stars and their basic stellar parameters has been extended to considerably larger 
data sets. The most significant contribution to the absolute magnitude Mv comes from the period P and from 
the first Fourier amplitude Ai, but there are statistically significant contributions also from additional higher 
order components, most importantly from A3 and in a lesser degree from the Fourier phase if 51. When different 
colors are combined in reddening-free quantities, we obtain basically period— luminosity— color relations. Due to 
the logTcff (B — V, log [Fe/H]) relation from stellar atmosphere models, we would expect some dependence also 
on Lp3i. Unfortunately, the data are still not extensive and accurate enough to decipher clearly the small effect 
of this Fourier phase. However, with the aid of more accurate multicolor data on field variables, we show that 
this Fourier phase should be present either in 1/ — / or in _B — V or in both. From the standard deviations of 
the various regressions, an upper limit can be obtained on the overall inhomogeneity of the reddening in the 



individual clusters. This yields oe(b- 
of (Jv>0.018 mag. 



V) 



:0.012 mag, which also implies an average minimum observational error 
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^ 1. Introduction 



In recent years we have conducted a series of studies aimed 
at deriving empirical relations between the Fourier param- 
eters of the light curves and the physical parameters of 
the RR Lyrae stars (Jurcsik & Kovacs 1996; Kovacs & 
Jurcsik 1996, 1997; hereafter JK96, KJ96, KJ97, respec- 
tively) . The method is purely empirical, and based on the 
assumption that the period and the shape of the light 
curve are directly correlated with the physical parameters 
(or quantities related to them), such as the absolute mag- 
nitude Mv, intrinsic color [B — V)q and metal abundance 
[Fe/H]. Therefore, it is hoped that once the formulae are 
derived for a representative sample, the relative values of 
the above physical parameters can be determined for any 
other RR Lyrae star, assuming that accurate light curve 
parameters in V color are available. 



Send offprint requests to: G. Kovacs 

* The National Optical Astronomy Observatory is operated 
by the Association of Universities for Research in Astronomy, 
Inc., under cooperative agreement with the National Science 
Fundation. 



The advantage of this method would be the utilization 
of very accurate observables (period, Fourier parameters) 
in calculating individual stellar parameters. If it were pos- 
sible to fit a large number of calibrating objects (in the 
present study, globular cluster variables) with the accu- 
racy of the observational noise, this would imply a precise 
calculation of the relative physical parameters from the 
derived formulae. In addition to the more practical appli- 
cations of relative distance and reddening determinations, 
the ultimate goal is to reach the level of accuracy at which 
a direct (star by star) comparison with the evolutionary 
calculations becomes possible. 



In JK96 we used Galactic field RR Lyrae stars for 
the derivation of the [Fe/H] formula. The expression was 
shown to give reliable abundances also for cluster stars, 
except perhaps at the low abundance end, where our for- 
mulae predict somewhat higher abundances. (However, we 
note that this result is based on the average cluster val- 
ues obtained from direct abundance analyses of giants and 
not of RR Lyrae stars. Recent spectroscopic investigations 
of Behr et al. (1999) indicate that substantial abundance 
differences might exist even within the horizontal branch.) 
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The relations for My and {B~V)o were obtained from 
globular cluster stars (see KJ96 and KJ97). If we assume 
that cluster and field RR Lyrae stars cover overlapping 
evolutionary stages and physical parameter space, then 
these latter formulae can also be applied to the field stars. 
The validity of this statement depends crucially on the 
completeness of the calibrating sample. Recent CCD pho- 
tometric studies of cluster RR Lyrae stars allow us to sat- 
isfy this condition more closely in the present investigation 
than in our former studies, and revisit the problem of ab- 
solute magnitude and color calibrations. Both the quality 
and the quantity of these new data lead to better deter- 
mination of the number of significant Fourier parameters 
entering in the formulae and, at the same time, to a more 
accurate calculation of the regression coefficients. 

2. The new empirical formulae 

The clusters and the corresponding number of variables 
together with the sources of the data are summarized 
in Table 1. This table contains all fundamental mode 
RR Lyrae (RRab) stars from the given sources, except 
those with: (a) obvious amplitude modulation (Blazhko 
effect), (b) clear outlier status due to e.g., cluster non- 
membership, (c) poor quality light curves, (d) short pe- 
riods and sinusoidal light curves (e.g., V73, 76, 169, 185, 
189 in uj Cen and V70 in M3). Fourier parameters of the 
V light curves and average colors used in this paper are 
given in Table 2.|^ Color indices and average V colors are 
calculated from the same source, except for M5. In the 
case of this cluster we average the V values used by KJ97 
and the ones published by Kaluzny et al. (2000) for all 
common stars in the two publications. In comparison with 
KJ96 and KJ97, we see that there is more than a factor 
of two increase in the number of stars in all colors (this is 
so even though in the present data sets we have omitted 
all photographic and less accurate CCD data used in our 
previous analyses). The large sample enables us to deal 
with separate data sets and obtain further information on 
the statistical significance of the various regressions. 

For deriving the relations between the Fourier param- 
eters, absolute magnitudes and colors, we follow the same 
methodology as in our former studies. For example, in the 
case of My , the best regression is obtained from the min- 
imization of the following expression 

j = l i=l k=l 

where N* = N — — Np, N is the total number of 
stars, Nf, is the number of clusters, Np is the number of 
the Fourier parameters, Nj is the number of stars in the 
j-th cluster, Vj{i) is the observed average magnitude of 
the z-th star in the j-th cluster. The Fourier parameters 
(including the period) are denoted by J-k{i)- The least- 
squares minimization leads to the determination of the 

^ Table 2 is available only in electronic form at CDS (ftp 
130.79.128.5) 



Table 1. Summary of the globular cluster data for RRab 
stars 



Cluster 


Nb 


Nv 


Ni 


Source 


M2 


13 


13 




LC99 


M3 




28 




K98 


M4 


4 


5 


3 


KJ97 


M5 


12 


43 


24 


KOO, KJ97 


M9 




4 




CS99 


M55 


4 


4 




099 


M68 


5 


5 


5 


W94 


M92 


5 


4 


3 


KJ97 


M107 




7 




CS97 


NGC1851 


11 


11 


11 


W98 


NGC5466 


7 


7 




C99 


NGC6362 


12 


12 


10 


W 


NGC6981 


20 


20 


18 


w 


IC4499 


49 


49 


41 


WN96 


Rup. 106 


12 


12 




K95a 


UJ Cen 




44 




K97 


Sculptor 




90 




K95b 


NGC 1466 


8 


8 




W92b 


NGC 1841 


9 


9 




W90 


Reticulum 


8 


8 




W92a 


Total : 


179 


383 


115 





References: CS97: Clement & Shelton (1997); CS99: Clement 
& Shelton (1999); C99: Corwin et al. (1999); K95a: Kaluzny et 
al. (1995a); K95b: Kaluzny et al. (1995b); K97: Kaluzny et al. 
(1997); K98: Kaluzny et al. (1998); KOO: Kaluzny et al. (2000); 
KJ97: Kovacs & Jurcsik (1997); LC99: Lee & Carney (1999); 
099: Olech et al. (1999); W90: Walker (1990); W92a: Walker 
(1992a); W92b: Walker (1992b); W94: Walker (1994); WN96: 
Walker & Nemec (1996); W98: Walker (1998); W: Walker (un- 
published) ; 

relative reddened distance moduli dj and the regression 
coefficients . The zero point of the magnitude scale mo 
is absorbed in the distance moduli and should be deter- 
mined from some direct distance calibration (e.g., paral- 
laxes, Baade-Wesselink analyses, double-mode variables, 
etc., however, please note that the absolute magnitude 
scale of the RR Lyrae stars is still an unsettled issue — 
e.g., Stanek et al. 2000 vs. Feast 1999 and Kovacs 2000). 

For each fixed Np we check all possible combinations 
from the period P and from the first five Fourier ampli- 
tudes and phases {Ai, A^; ip2i, y^si — see Simon & 
Teays (1982) and JK96 for the definition of these quanti- 
ties). After finding the best fit, we repeat the procedure 
with another Np (< 8). The optimum fit is obtained at 
the lowest parameter number for which the fitting accu- 
racy T) starts to become constant. Throughout this paper 
we use linear combinations of the Fourier parameters. We 
have not found any clear sign of preference for using non- 
linear combinations either for the fitting or for the fitted 
quantities. 

Three observables, related to static (or more accu- 
rately, to pulsation-cycle-averaged) stellar parameters, are 
studied. The first one is the intensity averaged V magni- 
tude. The advantage of this quantity is that it is directly 
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related to the absolute magnitude My of the static star 
(cf. Bono, Caputo & Stellingwerf 1995). The disadvan- 
tage of the V magnitude is that it is strongly affected 
by inhomogeneous cluster reddening which may be im- 
portant in some clusters. From this point of view, a bet- 
ter approach is fitting reddening- free quantities, such as 
W = V- Rv{B -V) imd X - Ri{V - I), where Ry 
and Ri are the extinction ratios (e.g., Cardelli, Clayton & 
Mathis 1989; Liu & Janes 1990). Here we set Ry = 3.1, 
Ri = 2.5 (the latter value corresponds to the Kron- 
Cousins system in / color). We use magnitude averages 
both in W and in X, because there does not seem to be a 
strong preference toward other averages when a compari- 
son is made between the static and average color indices 
in the nonlinear pulsation models (Bono ct al. 1995). The 
observed and intrinsic quantities are related through the 
true distance modulus 

W = Wo + d , (2) 

and similar expression holds also for X. Our task is to de- 
rive relations between the Fourier parameters and the in- 
trinsic quantities. For the calculation of the distance mod- 
ulus in the case of other values of Ry and Ri than the ones 
used in this paper, the expressions derived here for Wq and 
Xq can still be used but the distance modulus should be 
modified according to KJ97. 

The price we have to pay for eliminating reddening 
in W and X is the amplification of the observational 
noise due to the large values of Ry and Ri. If the red- 
dening is sufficiently small, it is obvious that using W 
and X is not the best idea, because of the unnecessar- 
ily amplified observational noise. Therefore, for any noise 
in the reddening and in the B and V magnitudes, there 
is an optimum trade-off between filtering out reddening 
fluctuations and keeping the effect of observational noise 
at minimum. Indeed, the expression for the variance of 
W = V — a{B — V) shows this property 

air = {Ry - afa% + o?(j% + + afa^ - 

- 2a{l + a)Kav(jB ■ (3) 

Here as, (Jb and ay are the standard deviations of the 
reddening, B and V magnitudes, respectively, and K is the 
correlation coefficient of V and B. Although tests with 
the observed data and the study of the above equation 
indicate that aw 1.5 - 2.0 would be more preferable than 
a = Ry, the gain is only about 10% in the increase of 
the fitting accuracy. This improvement does not seem to 
justify the giving up of the more transparent approach 
based on standard extinction ratios. 

A general problem in searching for the best-fitting lin- 
ear combination in the form of Eq. (1) is that the Fourier 
parameters J^k{i) are not independent. For the period and 
the total amplitude this correlation has been known for 
a long time (Preston 1959). As is shown in Fig. 1, there 
are approximate correlations also with the other Fourier 
parameters. Note that the scatter in these diagrams has 
a physical origin, because the accuracy of the low-order 



Table 3. Data sets for the My fit 



A 


B 


C 


N = 191 


N = 192 


N =383 


Nc = 15 


Nc = 5 


Nc = 20 


M2, M3, M4, M9, 


M5, uj Cen, 


A & B 


M55, M68, M92, M107, 


N5466, Sculptor, 




Rup 106, N1851, 


Reticulum 




N6362, N6981, IC4499, 






N1466, N1841 







Fourier parameters is very high (better than 1% if only 
the observational errors are considered, but may be less 
accurate than this in the case of hidden Blazhko variables 

— see also Sect. 2.1 for the discussion of the total errors in 
the Fourier parameters). One of the reasons for the scatter 
in Fig. 1 is that there are metallicity differences between 
the clusters. Assuming constant [Fe/H] for a given clus- 
ter, we would expect straight lines shifted vertically in 
the P (fsi diagram (see JK96). Because of the rela- 
tively narrow range, overlapping distributions of [Fe/H] 
in different clusters and chemical inhomogeneity in sev- 
eral clusters, this ridge structure is not seen in the figure. 
Other phases show similar progression, due to their tight 
interrelations with (^31 (JK96, see Kovacs & Kanbur 1998 
for an update of these relations). The scatter in the am- 
plitudes is only partially due to metallicity effects (they 
yield much weaker correlation with [Fe/H] than the phases 

— see JK96). Because of the significant amplitude depen- 
dence of the luminosity and temperature, the rest of the 
scatter is attributed to the variation of these quantities. 

The problem of the internal correlation of the Fourier 
parameters can be circumvented by constructing an or- 
thogonal set with the aid of principal component anal- 
ysis (sec Kanbur, Mariani & lono 2000). Although this 
approach certainly results in a better posed numerical 
scheme, both methods should end up with the same con- 
clusion, because both of them use linear combinations of 
the same set of Fourier parameters. More sophisticated 
methods, such as the one mentioned, could be instrumen- 
tal in the future studies of the finer (high frequency) de- 
tails of the light curves, when the Fourier decomposition 
becomes a non-useful concept. 

2.1. The absolute magnitude My 

Because we have a large number of stars observed in V 
color, we are in the position to test the significance of the 
derived formulae on two independent data sets. In Table 3 
we show the distribution of the clusters in the two sets. 
The clusters are sorted in the two groups by following 
the guideline of yielding about the same number of stars 
and including good quality large samples in both sets. In 
searching for the best fitting formula, we successively leave 
out those stars which can be regarded as outliers. This 
could be a somewhat delicate procedure, especially in the 
case of small sample size or noisy data, such as the sam- 
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Fig. 1. Fourier parameters of cluster RRab stars 



pies used in our former investigations. In this work we are 
more conservative than in our previous studies and confine 
ourselves to secure outliers which exceed the « 3cr limit. 
The list of these stars is given in Table 4. We see that 
there are 17 variables which do not conform to the rela- 
tion required by the rest of the sample. This means a 4% 
chance that the derived formula yields incorrect results in 
a randomly chosen sample of RRab stars (assuming that 
the sample does not contain 'obvious' Blazhko variables 
or stars with blended components). The chance of error 
may be even smaller, because some of the discarded stars 
might be Blazhko variables which would have been omit- 
ted if we had longer photometry available. Some of the 
stars listed in Table 4 show reasonably clear signs of ob- 
servational or other defects. For example, V20 in N1851 
and V16 in M2 have excessive scatter in their light curves. 
Blending or unresolved close components might be sus- 
pected in V118 of uj Cen and in V01446 and V02575 of 
Sculptor. Additional, more detailed examination of these 
and other stars are needed to clarify the reasons for their 
outlier status. 

Fig. 2 shows the significance of the derived formulae. 
The standard deviations at Np — are equal to the av- 
erages of the standard deviations of the observed mag- 
nitudes in each cluster. We see that the decrease of the 
unbiased estimates of the standard deviations for param- 
eter numbers greater than 2, is very small. Therefore, the 
present data seem to suggest the existence of a relation 
which contains only 2 parameters. This is in confiict with 
our previous result, where we argued for the validity of a 
3 parameter formula. Before we draw a conclusion, it is 



Table 4. Outliers in the various regressions 

Regression Outliers 

Mv M2: V4, 16, M5: V27, 59, M107: VIO, 

id Cen: V109, 118, 150, 154, N1851: V20, 
Sculptor: V01446, 02575, 06032, 06034, 
IC4499: V76, N1466: Wel4, N1841: V19 

Wo M2: V16, M5: V59, 86, M92: V4, 

N1851: V20, IC4499: V167, Rup. 106: VIO 

Xo M5: V13, 27, 59, 83, N1851: V20, 

IC4499: V34 



o 
o 

+ 




Fig. 2. Standard deviation of the residuals of the My fit 
as a function of the number of the regression parameters. 
Results displayed for data sets A and B are vertically 
shifted for better visibility {const. = 0.048, 0.015 and 0.0 
for sets A, B and C, respectively 

necessary to test if the smaller change of V in the case of 
the high parameter fits is statistically significant. Because 
of the complexity of the fitting procedure (e.g., the use of 
noisy Fourier parameters), it is not possible to apply stan- 
dard methods of analysis of variances (e.g., Sachs 1982). 
Instead, direct numerical tests are performed. More specif- 
ically, in testing data set C, the procedure below is fol- 
lowed. 

For each parameter number Np, synthetic data are gen- 
erated by using the formula obtained from the given data 
set. Then, Gaussian noise is added both to these data and 
to all Fourier parameters i.e., the test data have the fol- 
lowing form: 



E 

k=l 



mk[Tk{i) + ?7fc(i)] 



(4) 



In these tests the role of the exact values of the relative dis- 
tance moduli dj is not essential. Therefore, we take simply 
the average V magnitudes in each cluster as very rough 
substitutes for the exact values. The noise components are 
chosen in such a way as to yield the observed standard 
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Table 5. Statistical test of the significance of the various 
My regressions on data set C (see text for details) 





Pi{%) 


_.max fO/ \ 


^obs 






S 


1 


100 


7 


1541 


19 


18 


85 


2 


100 


6 


139 


15 


13 


9 


3 


67 


27 


47 


11 


9 


4 


4 


62 


27 


8 


10 


9 





5 


36 


55 


8 


8 


8 






Note: Ao6s, Atest and ua are in 10 ® mag^ 



deviation at the given Np. Because the observed standard 
deviation is a result of the errors of the Fourier parame- 
ters and of the average magnitudes, is slightly smaller 
than the observed standard deviation. Unfortunately, the 
accurate values of the errors of the Fourier parameters are 
not known. Although the formal statistical errors can be 
computed, because of the larger dispersion of the interre- 
lations, there might be additional sources of errors (e.g., 
Blazhko effect, background contamination). From the up- 
dated relations of Kovacs & Kanbur (1998) we think that 
the values of a{Ai) = 0.01 and cr((^ji) = 0.03 are reason- 
able estimates for the overall total errors of the low-order 
Fourier parameters. These values are used at all Np in the 
following simulations. The adjusted standard deviations 
of are 0.057, 0.041, 0.039, 0.038 and 0.037 for Np = 1, 
2, 3, 4 and 5, respectively. We note that these standard 
deviations are up to 0.002 smaller than the observed ones 
shown in Fig. 2. 

For each realization of the data generated above, stan- 
dard parameter searches are performed, similar to the ob- 
served data. By repeating the analysis for many realiza- 
tions, information is obtained on the significance of the 
regression at parameter number Np. 

The following statistical parameters are calculated: 

— Pi: probability of the identification of the parameter 
set used in the generation of the synthetic data. 

— cr™"^: maximum of the relative standard devi- 
ations of the regression coefficients (cr™"^ = 
max{a{mk)/\mk\}). 

— A = I>2(7Vp) - I>2(A/p + 1) and its standard devia- 
tion (TA- 

— S = (Ao6s — Ate5t)/<JA, the significance of the A^p -|- 1 
versus the Np parameter fits. 

The results of the above test on data set C are shown 
in Table 5. We see that the two parameter regressions 
are undoubtedly very significant in all test quantities. 
Furthermore, a little closer examination of the three and 
four parameter regressions leads to the conclusion that 
these higher parameter formulae are also significant (al- 
though the four parameter one might be close to the detec- 
tion limit). The present data render any further increase 
of the parameter number to be statistically insignificant. 
This result is basically in line with the conclusion of KJ96, 
although here we get A3 and (^51 for the higher param- 
eters, whereas in KJ96 we obtained <^3i. This difference 
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Table 6. Formulae for My obtained from various data 

sets 



Set Mv + const = V 



A 


-1.682 log P 


- 0.736yli 




0.0365 


B 


-1.897 log P 


- 0.855Ai 




0.0460 


C 


-1.820 log P 


- O.8O5A1 




0.0416 


c 


-1.876 log P 


- 1. 158^1 


+ O.82IA3 


0.0399 


c 


-1.963 log P 


- 1.124.4i 


+ 0.830.4.3 










+ ().Oii^.-,i 


0.0393 



Notes: — phase ipsi must be in the [3.5,5.5] interval 

— outliers have been omitted according to Table 4 



is attributed to the smaller data set used in KJ96 and to 
the interrelations between the Fourier parameters, already 
mentioned at the beginning of this section. We return to 
the compatibility of this work and that of KJ96 in Sect. 4. 

In column 5 the i^test values are positive, although we 
would expect zero, since at any Np the higher parame- 
ter regressions fit only the noise, which should give zero 
Afesi on the average, because this quantity is calculated 
from the unbiased estimates of the standard deviations 
(see Eq. (1)). The reason why we get small positive values 
is that in the Np + 1 parameter search, those parameter 
combinations are selected which yield the smallest disper- 
sion from the ones which employ Np + 1 parameters from 
the first 10 Fourier components. Therefore, there is a bet- 
ter chance to get some decrease in the dispersion then if 
we had only one additional component to use in the Np-\-l 
parameter fit. 

We note that if we did not use noise in the Fourier pa- 
rameters, we would reach the same conclusion about the 
parameter number, although with an even higher signifi- 
cance. 

Formulae obtained for My with various parameter 
numbers from data sets A, B and C are given in Table 6. 
The 10% difference between the coefficients of log P of the 
formulae derived on sets A and B can be accounted for 
by the errors in the coefficients (the Icr relative errors of 
the regression coefficients are f« 4% for both data sets). 
The two formulae yield also very similar My values. A 
comparison made on a sample of w 500 stars compiled 
from Galactic field and cluster variables, gave 0.011 mag 
for the standard deviation of the differences between the 
My values computed from these formulae. Considering 
the total range of My (which is « 0.6 mag for this large 
data set and « 0.4 mag for the clusters), this standard de- 
viation corresponds to compatibility at the 2% level. For 
data set C we see changes in the coefficients with the vari- 
ation of the parameter number. This is primarily due to 
the interrelations among the Fourier parameters. As is ex- 
pected, the differences between the various My estimates 
obtained from the high- and low-parameter formulae are 
rather small. By using the above large data set, we get 
a = 0.012 and 0.008, if we compare the two and three, and 
three and four parameter formulae, respectively. Since — 
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-.1 .1 .2 .3 .4 .5 



Vo(P, Af, A3) + const. 

Fig. 3. Observed V magnitudes versus absolute magni- 
tudes (Vo, calculated from Eq. (5)). Lines show fiducial 
values obtained from Vq after shifts with optimum dis- 
tance moduli. Inset displays the same quantities after ad- 
justment of the V magnitudes by the relative distance 
moduli 



as we have already seen — the three and possibly even 
the four parameter formulae are significant, this compari- 
son gives some indication on the accuracy to be gained by 
the application of the high parameter formulae (however, 
this gain is partially counterbalanced by the statistical er- 
rors of My — see Appendix A). This level of change in 
My is not very important in simple applications (e.g., in 
distance modulus calculations), but could have interest- 
ing consequences when empirical evolutionary tracks are 
attempted to be calculated. 

Fig. 3 displays the correlation between the observed 
and calculated magnitudes. (It is noted that M107 and M4 
are not plotted. For M107 the reddened distance modu- 
lus is equal to that of M3. The cluster M4 is omitted, 
because plotting would result in a large crowding in the 
figure due to its high observed brightness.) The calculated 
magnitudes are obtained from the three parameter fit to 
the magnitude averaged values. This yields the following 
formula for data set C 

yo = -l-8801ogP- 0.971^1 + O.9OM3 + const. . (5) 



The slight difference between this and the corresponding 
formula in Table 6 is due to the different ways of averag- 
ing of the observed magnitudes. The above formula fits the 
data with a — 0.040 mag. The nice correlation with the 
variables of the individual clusters is very comforting. This 
shows that with a good chance, no important physical de- 
pendence in Vo (or in My) are missed in the above deriva- 
tion. The method could have been a failure if My were a 
single function of [Fe/H] and, at the same time, the clus- 
ters were chemically homogeneous. This situation would 
lead to 'truly horizontal branches', which would mean the 
absorption of all physical parameters into the (incorrect) 
distance moduli and the lost of any correlation with the 
Fourier parameters. Both chemical inhomogeneities and 
complex parameter dependence of My prevent this sit- 
uation from happening. Since most of the clusters have 
reasonably small [Fe/H] dispersion, we can conclude that 
the empirical [Fe/H] — My relation is rather loose, as has 
already been shown in KJ96. 

2.2. The PLC relation for B, V 

We use the data sets given in Table 7 to find the optimum 
formula for the reddening- free quantity W. By following 
the same procedure as in the My fit, outliers are omit- 
ted iteratively. The list of the 7 stars which were left out 
from the final samples is given in Table 4. On the basis of 
the small number of outliers we find similar 'applicability 
ratio' of the derived formulae as in the case of the My 
fit. The variation of the fitting accuracy as a function of 
the parameter number is shown in Fig. 4. The single pa- 
rameter nature of W seems to be very strongly indicated. 
This result is at variance with the conclusion of KJ97, 
who argued for the significance of a 3 parameter relation. 
The reason of this contradiction is that KJ97 used con- 
siderably smaller sample with larger observational noise 
(e.g., photographic data for N3201 and M107). This re- 
sulted in a somewhat biased selection of the outliers, and 
a concomitant overestimation of the parameter number. 
Nevertheless, as we shall see in Sect. 4, the derived formu- 
lae give similar results. 

Table 8 shows that data sets a and b yield almost the 
same single parameter formulae. In the last row we also 
display the two parameter formula which might be sus- 
pected as the highest parameter relation allowed by the 
present data set for Wq. Although the absolute change in 
2? is similar to the one observed between the three and 
four parameter formulae for My, here the statistical sig- 
nificance is even lower, because of the more than a factor 
of two lower number of data. Indeed, repeating the same 
type of numerical test as in the case of My , we get a sig- 
nificance ratio of S — 1.7. Furthermore, for the occurrence 
rate Pi of the two parameter formula and for the largest 
parameter error cr™°^ we obtain 55% and 39%, respec- 
tively. All these results leave only a rather slim chance for 
the statistical significance of the two parameter formula 
for Wq. In addition, if the predicted Wq values are com- 
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Table 7. Data sets for the Wq fit 
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a 


b 


c 


N =87 


N =92 


N = 179 


iVc = 10 


Nc = 5 


iVc = 15 


M4, M5, M55, M68, 


M2, N1851, 1C4499, 


a & b 


M92, Rup. 106, N6981, 


N6362, N5466 




Retic, N1466, N1841, 







Table 8. Formulae for Wq obtained from various data sets 



Set 


Wo + const = 


■D 


a 


-2.477 logP 


0.0434 


b 


-2.462 logP 


0.0380 


c 


-2.467 log P 


0.0405 


c 


-2.204 logP + 0.182Ai 


0.0397 



Note: — outliers have been omitted according to Table 4 




1 2 3 4 5 6 7 



Fig. 4. Standard deviation of the residuals of the W fit 
as a function of the number of the regression parame- 
ters. Results displayed for data sets a and b are vertically 
shifted for better visibility {const. = 0.036, 0.018 and 0.0 
for sets a, b and c, respectively 

pared on the large data set mentioned above, they agree 
with a — 0.014 mag standard deviation. This, with the 
total ranges of 0.8 and 0.7 mag for Wq on the field & 
cluster, and cluster variables, respectively, corresponds to 
compatibility at the 2% level. Therefore, in the following 
we consider only the single parameter formula as a statis- 
tically significant relation. 

The strong correlation of the period with W for the 
individual cluster variables is shown in Fig. 5. As for the 
Aly regression, we emphasize the importance of this cor- 
relation from the point of view of the separation of the 
dependence of W on the period, from that of on the dis- 
tance modulus. 

The single parameter dependence of Wq establishes a 
period-luminosity- color (PLC) relation for RRab stars. 



21 
20 
19 
18 
17 




16 



15 
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11 



log P [day] 

Fig. 5. Observed reddening-free quantity W — V ~ 
3.1{B — V) versus period. Lines show fiducial values ob- 
tained from Wq (see the single parameter formula for data 
set c in Table 8) after shifts with optimum distance mod- 
uli. Inset displays Wq versus W — d, where d is the relative 
distance modulus 



The existence of this type of relation is very basic for 
Cepheids (e.g., Udalski et al. 1999). On the other hand, for 
RR Lyrae stars, to the best of our knowledge, this is the 
first time that the existence of this relation is exhibited 
and an accurate formula is given. Other works devoted to 
the study of the relations between the main physical and 
observed parameters of RR Lyrae stars have also revealed 
certain aspects of the Aly and Wq relations derived in this 
paper. For example, Nemec, Nemec & Lutz (1994) found 
My - \ogP - [Fe/H] relations, whereas Castellani & De 
Santis (1994) claimed the existence of a logTcff - logP - 
Ab (blue amplitude) relation. Sandage (1981) derived a 
P - L - A relation by combining theoretical and empiri- 
cal data. The results given in this paper are distinct from 
the above ones, because our formulae have been derived 
by using large data sets, which arc based on the best cur- 
rently available observations. In addition, our approach is 
completely empirical, without resorting to theoretical as- 
sumptions and few parameter (e.g., period, B amplitude) 
searches. As a result, the formulae are more accurate. For 
example, the relative accuracy of the coefficient of log P 
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in the expression for Wo derived from data set c is 2.4% 
(for a complete error formula see Appendix A). 

From the formula of Table 8 and Eq. (5), one can derive 
a relation for the intrinsic color index {B — V)o. When the 
zero point is fitted with the one established for {B — V)o 
by KJ97, we get 

{B - V)o = 0.189 logP - 0.313^1 + 0.293^3 + 0.460 . (6) 

A more conservative expression is obtained if the two pa- 
rameter counterpart of Eq. (5) in conjunction with the 
above single parameter formula for Wq is applied 



iB-V)o = 0.2091ogP-0.187v4i + 0.453 



(7) 



This second formula is given to aid those applications in 
which only limited information is available on the light 
curves and when the goal is to compute average quanti- 
ties on large samples or when one needs to combine less 
accurate data with the above expressions (e.g., reddening 
estimations). In any case, Eqs. (6) and (7) yield very simi- 
lar results. By repeating the same compatibility test as for 
My, we obtain a — 0.004 mag for the standard deviation 
of the differences of the {B — V)o values obtained from 
these two equations. With the total {B — V)o range of 
0.11 mag for the RRab stars, this result implies compati- 
bility at the 4% level. A similar conclusion can be drawn 
if a comparison is made between the present formulae and 
the one given in KJ97 (see Sect. 4 for details). 

2.3. The PLC relation for V, I 

In KJ97 we made an attempt to derive a formula also for 
the other reddening- free quantity X, which contains V and 
/. However, due to the small amount of direct cluster data 
available then, we combined those data with the ones for 
field stars, through the application of the Wq and {B — V)o 
formulae. This is certainly not an independent derivation 
of the formula for Xq. Here we are in the position to give 
a direct estimate on Xq, based solely on cluster variables. 

Because the available data are still not too extensive, 
we cannot make tests on independent data sets as we did 
in the case of Wq. By proceeding exactly in the same way 
as in the case of the W fit, we obtain the result shown 
in Fig. 6 (for the list of the outliers, see Table 4). It is 
clear that the single parameter dependence of Xq is very 
strongly suggested. For the 109 variables we get 



Xq = -2. 513 log P + const. 



(8) 



This expression fits the data with cr = 0.037 mag. The 
error of the coefficient of log P is only 3% (see Appendix A 
for the complete error formula). To exhibit the tightness 
of the correlation of Eq. (8) with the observations, Fig. 7 
displays the fit to the variables of the individual clusters. 

Due to the single parameter dependence of Xq, the 
present result is in conflict with the one given in KJ97. We 
recall the same facts as in the case of Wq to explain this 
contradiction (see also next section for further discussion 
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Fig. 6. Standard deviation of the residuals of the X fit as 
a function of the number of the regression parameters 



18 



17 



16 



15 



14 



13 



12 



11 



10 




IC4499 



N1851 



N6362 




N6981 



M92 



log P [day] 

Fig. 7. Observed reddening-free quantity X = V~2.5{V~ 
I) versus period. Lines show fiducial values obtained from 
Xq (Eq. (8)) after shifts with optimum distance moduli. 
Inset displays Xq versus X — d, where d is the relative 
distance modulus 



of the question related to the number of parameters in Xq 
and Wo). 

In the same way as in the case of (-6—^)0, the following 
expressions can be derived for the intrinsic color index 

(y - /)o = 0.253 log P - 0.388^1 + 0.364^3 + 0.648 , (9) 
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or, by using the two parameter formula for Vq 
(l/-7)o = 0.2781ogP- 0.232^1 + 0.640 . 



(10) 



Assuming standard extinction ratios (sec Sect. 2), the 
zero points in these formulae are consistent with those of 
{B — V)o in the sense that they yield the same reddening 
with (TAEiB-v) = 0.023 for the clusters with simultaneous 
B, V and / observations {AE{B — V) denotes the differ- 
ence between the reddenings obtained from the B — V and 
V — I indices). It is suspected that at least part of this 
large scatter is due to zero point errors of the various color 
indices in the individual clusters. Color indices calculated 
from Eqs. (9) and (10) agree with a = 0.005 mag. 

3. The problem of T^s consistency 

The derived expressions for the dereddened colors can be 
utilized in calculating effective temperatures. By using the 
stellar atmosphere models of Castelli, Gratton & Kurucz 
(1997), simple Unear formulae can be obtained for logTeff 
as a function of the color index, gravity log 5 and metal 
abundance [M/H]. For further reference, here we repeat 
these formula from Kovacs & Walker (1999) 

logTeff = 3.8840- 0.3219(5- y)o + 0.01671ogc/ + 

+ 0.0070[M/H] , (11) 

logTeff = 3.9020 -0.2451(F-/)o + 0.0099 log 5 + 

- 0.0012[M/H] . (12) 

The range of validity of these formulae covers the pa- 
rameter regime of cluster RRab stars, in particular, 
—2.5 <[M/H]< —0.5. Therefore, the various temperature 
estimates are to be compared on the set of 383 cluster 
variables with V observations (see Table 1). It is men- 
tioned that while Eq. (12) fits the corresponding model 
values with (T(logTeff) = 0.001, for Eq. (11) this accu- 
racy is only 0.003. Better fitting formulae with {B — V)o 
can be obtained by choosing narrower parameter ranges. 
Nevertheless, these formulae are very similar to Eq. (11), 
and therefore, it is not surprising that tests made with 
them have led to the same conclusions as the ones to be 
discussed below. 

By using Eqs. (6) and (9), the [Fe/H] formula of JK96 
and the expression given for logg by Kovacs & Walker 
(1999), we can compute the temperatme for the two color 
indices and compare the two estimates. The result is 
shown in Fig. 8. A similar correlation is obtained if Eqs. (7) 
and (10) are used for the color indices. The standard devi- 
ation of the logTeff differences is 0.0026 mag. Considering 
the total range of « 0.04 mag, this size of dispersion does 
not seem to be small enough, although it is worth recalling 
the accuracy of the Tes values obtained from the B — V 
index (see above). More importantly, it is well known and 
is also clearly exhibited in Eq. (11) that Tes{B — V) has 
a strong [M/H] dependence. Since this quantity is a func- 
tion of we would expect some dependence on this 
Fourier phase to be observable at least in one of the color 
indices. Although the effect is small, the absence of 1^31 is 
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Fig. 8. Relative effective temperatures calculated from 
color indices B — V and V — I. Eqs. (6) and (9) are used 
on the full V data set for cluster variables (see Table 1). 
Filled circles denote variables with [Fe/II]> -0.5. The 45° 
line is shown for reference 



Or 




Qcalc. 

Fig. 9. Observed vs. calculated reddening-free quantity 
Q ^V -I - 1.24(5 - V). Cluster variables are plotted by 
open circles, field variables by filled circles. The 45° line 
is shown for reference 

somewhat puzzling. We suspect that noise, limited [Fe/H] 
range of the globular clusters and delicately distributed 
(fisi dependence between the two colors are responsible 
for this contradiction. In the following we show a deriva- 
tion of the relation between the two color indices, which 
might support the above conjecture. 

The idea is to avoid distance-dependent quantities, 
such as X, and utilize more accurate multicolor obser- 
vations. It is easy to see that the quantity Q = V — I — 
1.24(5 — V) is reddening-free (assuming standard extinc- 
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tion ratios — see Sect. 2), independent of the distance, and 
enables us to use the more accurate data on field RRab 
stars. By employing the sample of 25 variables from KJ97, 
with the omission of the outlier SS Leo, the following for- 
mula is obtained 

Q = 0.477 + 0.397 log P- 0.060^331 • (13) 

This expression fits the 24 variables with a = 0.011 mag. 
The single parameter regression contains the period and 
yields a = 0.016 mag. Because of the low number of data 
points, higher order regressions are not considered, but 
they all settle down at around cr = 0.0095 mag. The rela- 
tive errors of the regression coefficients in the above equa- 
tion arc 7% and 18%, for the coefficients of log P and ip^i, 
respectively. 

Adding cluster variables to this sample increases the 
noise considerably. Fig. 9 shows the regression for this 
large sample. (The variables of M4 arc omitted, because 
of their specific extinction ratios — see KJ97 and refer- 
ences therein.) In spite of the large scatter, this sample 
of 121 stars yields a very similar formula, basically with 
the same parameters: Q = 0.431 -1-0.409 log P-0.051(p3i. 
The significance of the two parameter formula is also very 
similar to the one obtained on the sample of field stars. 
The large scatter in this figure draws the attention to the 
importance of accurate photometry if more sophisticated 
problems, such as the one here, are studied. The majority 
of the outlying points in the upper left and middle right 
parts of the figure correspond to variables from IC4499 
and N1851, respectively. Reddening could not be the cause 
of this scatter, because IC4499 has a relatively large red- 
dening of Eb-v = 0.22, whereas N1851 has a low one 
of Eb-v = 0.03. The cause of the discrepancy apparent 
in these clusters is not known. Less accurate photometry 
with zero point errors, as well as peculiar extinction ra- 
tios are suspected to be the main possible agents in the 
systematic displacements of most of the variables of these 
clusters. For example, zero point errors of w 0.02 mag in 
the color indices might easily cause such a discrepancy. 

Fig. 10 shows the improvement of the compatibility 
of the temperatures derived from the different colors, if 
Eq. (13) is used to calculate V — I. Although some of 
the metal rich stars have become 'overcorrected', there is 
a definite improvement, indicating that the predicted <p3i 
dependence from the theoretical temperatures and the one 
obtained from the above, completely independent deriva- 
tion agree fairly well. However, it is also noted that be- 
cause {V — I)o is not calculated independently, the tem- 
peratures in this figure are highly correlated, which also 
decreases the scatter. We mention that metal rich stars 
are very sensitive to the contribution of (p^i in the color 
relation Q. For example, switching to the relation for Q 
obtained from the larger sample (see above), metal rich 
stars become even less discrepant. 

In closing this section we emphasize that the above 
test has been used to demonstrate only that there is a 
hope to derive consistent sets of equations for both color 
indices, but it is not claimed at this stage that we found 
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MogT,„(B-V) 

Fig. 10. Relative effective temperatures calculated from 
color indices B — V and V — I. Eqs. (6) and (13) are used 
to compute the B — V and V — I color indices, respectively. 
The full V data set for cluster variables (see Table 1) is 
plotted. Variables with [Fe/H]> —0.5 are plotted by filled 
circles. The 45° line is shown for reference 

the final formulae. The more definitive statement about 
this rather delicate issue may come only when additional, 
more accurate cluster data will become available. 

4. Distance moduli, reddenings, comparison witli 
Kjge and KJ97 

Average values and standard deviations of Eb-v and 
those of the various distance moduli are summarized in 
Table 9. Abundances — computed by the formula of JK96 
— are also listed for completeness. Except for the calcu- 
lation of [Fe/H], data sets without the outliers, as given 
in Table 4 were used. For [Fe/H] we employed the full set 
C of Table 3. The outliers (which were discrepant in re- 
spect of the average cluster [Fe/H] values) are listed in 
the note to Table 9. These variables have been omitted in 
the calculation of the averages and standard deviations. 
Reddening has always been calculated from the B and V 
colors and not from the less extensive (and probably also 
less accurate) I data. Because of the still debated issue 
of the luminosity of the RR Lyrac stars, we show only 
relative distances, all normalized to zero for IC4499. 

As we have already noted in our former papers, and is 
also clearly seen in Table 9, in some clusters the derived 
mctallicity values exhibit considerable scatter. It is worth 
noting that recent spectroscopic analysis of the horizontal 
branch stars in M13 by Behr et al. (1999) also revealed 
substantial inhomogeneities. It is seen that most of the 
clusters have metallicities around —1.5 and there are only 
a few with [Fe/H]>-1.0 or [Fe/H]<-1.9. This distribution 
of [Fe/H] signals some warning in respect of the unbiased 
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Table 9. Metallicities, reddenings and relative distance moduli 
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Notes: — dy denotes reddened distance moduli, calculated from color V 

— Rv =4.1 and Ri — 3.3 are used in the calculation of the distance moduli for M4 

— a{EB-v) = <^[{B — V)obs — {B — V)o], i.e., a{EB-v) contains both observational noise and contribution from 
reddening inhomogeneities 

— Zero points of the distance moduli arc sot to yield zero distance moduli for IC4499 

— Outliers with respect to the average cluster abundances and dispersions (derived [Fe/If] values are shown in paren- 
theses): M2: V708(-0.14); M3: V50(-1.91), V106(-0.76); M5: V6(-1.71), V38(-0.58), V56(-0.66), V963(0.60); 
M107: V12(-0.06); N1851: V8(-1.63), V20(-0.66); N6362: V361(-0.41); N6981: V297(-0.86); IC4499: 
V30(-0.32), V61(-0.73), V72(-0.78); uj Cen: V77(-1.21), V88(-1.16), V114(-0.80), V137(-2.06), V146(-1.07), 
V154(-1.15); Rup. 106: V5(-0.58) 



coverage of the possible evolutionary stages by the present 
sample. Specifically, one should exercise some caution in 
extrapolating the results of this paper to more metal abun- 
dant variables, such as the ones observed in the Galactic 
field. 

Reddening could also be inhomogeneous, with an ob- 
served overall standard deviation of ~ 0.014 mag. Because 
observational noise is also present in the data, the true dis- 
persion of Eb-v is smaller than this formal value. A more 
accurate upper limit on the overall reddening dispersion 
can be derived with the aid of some basic equations con- 
necting various noise properties. Denoting the standard 
deviations of the residuals of the Vb, Wo and {B — V)o 
regressions by (t{V), a{W) and (t{B — V) respectively, it 
is easy to see that the following relations hold for these 
quantities 

a\V) = al + Rlal , (14) 

- 2Rv{l + Rv)K* , (15) 
a^{B -V) = al + al + al- 2K* , (16) 

where K* = KuvO'b, K and the other quantities have the 
same meaning as in Eq. (3). Elimination of ay and ob 



from these equations results in the following expression 
for ge 

2 ^ (1 + 2i?y)a"(y) + Rln^B -V)- u'^jW) _ 
~ 2R^{1 + Rv) 



Rv{l + Rv) ■ ^ ^ 

Although the correlation coefficient K is not known, it 
can be safely assumed that it is positive (e.g., random 
background contamination influences both colors in a sim- 
ilar way). From Tables 6 and 8 we get a{V) = 0.040 and 
a{W) = 0.041. A direct fit to the observed {B - V) values 
yields a{B - V) = 0.014. These values and the K* > 
condition leads to the upper limit of 0.012 mag for cje- This 
constraint yields also a lower limit for ay - From Eq. (14) 
we get (Ty > 0.018 mag. (This latter quantity is very sen- 
sitive to the value of as- The value quoted was obtained 
from the formal value of as = 0.0115.) These limits are 
average values, deviations from them might occur in the 
individual clusters. 

In comparing the dw and dx distance moduli, we 
see that substantial differences might exist between them. 
Unfortunately, these true distance moduli derived from 
medium wavelength optical photometry are rather sensi- 
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Fig. 11. Comparison of the empirical relations derived in 

this paper and those given by KJ96 and KJ97. For the 
employed relations see the three parameter formula for 
data set C in Table 6 and Eq. (6). Vertical axes: this pa- 
per; horizontal axes: KJ96 & KJ97. The zero points are 
arbitrary. The 45° lines are plotted for reference 

tive to zero point errors. Indeed, if the zero points are 
changed by AV, A_B and A/, assuming standard extinc- 
tion ratios, the difference between the derived distance 
moduli from W and X becomes dw — dx = 5.5AV — 
S.lAi? — 2.5A/. With the commonly admitted zero point 
errors of 0.01 mag, it is very easy to accumulate a differ- 
ence of 0.1 mag between the two types of distance moduli. 

As we have already noted at the beginning, the present 
data sets are much more extensive then the ones used in 
our previous studies. Therefore, it is a matter of interest 
to check the difference between the old and new formulae. 
For the following reasons, in this comparison we restrict 
ourselves only to My and {B — V)o. For the V — I index 
the data are still limited and the result of the compari- 
son depends on whether we rely solely on the fit obtained 
from cluster variables, or incorporate the data available 
on field stars (as it was done in KJ97). Therefore, it is not 
possible to make a consistent comparison. For {V — K)o, 
there are only field data available with sufficient accuracy. 
In addition, because of their limited number, the result 
depends on the selection of the outliers. Therefore, we do 
not deal with this color either. 

The comparison for My and (B — V)o is made on the 
same large data set of 500 stars which was already em- 
ployed for similar tests in Sect. 2. The results are dis- 
played in Fig. 11. We can conclude that both quantities 
were reasonably well approximated already in KJ96 and 
KJ97. 

5. Discussion and conclusions 

In an effort to utilize the steadily increasing number of 
light curves accumulated during the various photometric 
surveys and CCD works on individual clusters, we revisited 
the problem of absolute magnitude and intrinsic color cal- 
ibrations in terms of the Fourier parameters of the light 
curves in V color. The amount of data used in this paper is 
more than a factor of two larger than the one used in our 
similar studies conducted a few years ago. The increase of 



the data base is beneficial for the following reasons: (i) the 
larger data sets enable us to derive statistically more sig- 
nificant relations; (ii) the number of parameters entering 
in the finally adopted formulae does not depend on the 
omission of the outliers, whose treatment could be a del- 
icate problem in the case of small sample size; (iii) with 
the larger number of variables, it is hoped that there is 
a better representation of the various evolutionary stages 
and also physical parameters, most importantly luminos- 
ity levels and abundances. It is this latter point which is 
perhaps the most crucial in the applicability of the derived 
formulae. In the present sample most of the data concen- 
trate around [Fe/H]Ri —1.5 and only a few clusters occupy 
the extremes close to [Fe/II]= —1.0 and —2.0. Therefore, 
some degree of caution is necessary when extrapolating 
the present results to variables with different metallici- 
ties, and especially beyond the high metallicity end (e.g., 
metal rich Galactic field variables). 

The effective temperatures calculated from the expres- 
sions for B — V and V — I agree with (T{A\ogTcs) = 0.003. 
The lack of higher level compatibility is due to the absence 
of phase dependence in our formulae which is required by 
the explicit [Fe/H] dependence in Tes, predicted by stel- 
lar atmosphere models and by our formulae derived earlier 
for [Fe/H]. The test presented in Sect. 3 suggests that the 
limited [Fe/H] coverage of the cluster data together with 
the still sufficiently high observational noise in the cluster 
data, are probably the main sources of the lack of better 
agreement between the various temperature estimates. 

For the absolute magnitude My, a highly significant 
three parameter formula has been derived with period P 
and Fourier amplitudes Ai and A3. There is also an indi- 
cation for the presence of the Fourier phase (p^i, although 
with a much lower significance. In certain applications 
(e.g., in computing averages on large representative sam- 
ples) there are very little differences between the high and 
low parameter formulae. In the case of less accurate data, 
one can use even the two parameter formula as it is given 
in Table 6. 

For the reddening-free quantities W = V — Rv{B — V) 
and X = V — Ri{V — I) we obtained basically single 
parameter formulae. Although these PLC relations are 
rather tight and well defined, it is worth remembering that 
this is not entirely consistent with the requirement posed 
by stellar atmosphere models (see above). 

Although the accurate value of the observational noise 
is not known, from the standard deviations of the residuals 
of the various regressions we may have estimates both on 
the observational noise and also on the dispersion of the 
reddening. For the latter quantity an upper limit was ob- 
tained, with as < 0.012 mag. This value puts a lower 
limit on the average observational noise, giving ay > 
0.018 mag. From this result and the fit of the reddening- 
and distance-free quantity Q = V — I — 1.24:{B — V) we can 
conclude that although the presently available CCD pho- 
tometric data on cluster variables allow us to derive more 
accurate empirical relations than ever before, the average 
colors (mostly due to crowded field effects) may still be not 
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accurate enough to address more sophisticated questions, 
such as the higher level compatibility of the temperature 
estimates obtained from various colors. 

From the point of view of applicability, it is worth men- 
tioning that the estimated probability that the present 
formulae give incorrect (i.e., « Scr discrepant) result is 
less than 5% (assuming that they are employed on RRab 
stars with stable light curves, which cover approximately 
the same [Fe/H] range as the ones in our sample). The 
formal statistical accuracy of the individual estimates are 
better than cr(My) w 0.01 and a{Wo) ~ cr(Xo) w 0.005. 
Therefore, these formulae are suitable (under the condi- 
tion just mentioned) for the accurate calculation of rel- 
ative absolute magnitudes, distance moduli and redden- 
ings. When combined with stellar atmosphere models, 
the formulae are very useful in mapping the instability 
strip (Jurcsik 1998) and comparing with the predictions of 
pulsation models (Kollath, Buchler & Fcuchtinger 2000). 
Further improvement of the formulae requires more accu- 
rate cluster data together with a better sampling at the 
low and high metallicity limits. Considering the steady 
progress in cluster photometry, we think that these goals 
are reachable in the very near future. 
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Appendix A: 

In the following we give a summary of the statistical error 
formulae for My, Wq and Xq. Starting with the two pa- 
rameter expression for My (see the 3rd row of Table 6), 
assmning that the period is error-free, the variance can be 
written in the following form 

3 

= 0.805Vi^ + 5^ ^'^P'P^ ' ("^-1) 

where a = 0.0416, pi = 1, p2 = logP + 0.22488, ps = 
Ai — 0.32161, and the symmetric correlation matrix Kij 
is given in Table A.l. We see that the Ku elements arc 
zero (lower than 10~^), because the corresponding sample 
averages have been subtracted from the parameters. In the 
same way as above, for the three parameter expression of 
My (see the 4th row of Table 6), we get 

4 

al,^ = 1.158Vi^ + 0.821^^3 + ^^PiPi ' (^-2) 

i,i = l 

where a = 0.0399, pi — P3 are the same as above, p4 = 
As — 0.09926. Table A. 2 displays the correlation matrix. 
Because Wq and Xq depend only on the period, their error 
formulae arc simpler 

= 0.04052[0.00581 + 2. 113(logP + 0. 22395)2], (A.3) 
0-1^ = 0.03662[0.00917 -F4.311(logP-F 0.23329)2]. (A.4) 



Table A.l. Correlation coefficients K^j in the error for- 
mula Eq. (A.l) 
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j 
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Kij 
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1 


0.00273 
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1.63537 
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2 


0.00000 
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0.74685 


1 


3 


0.00000 
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3 


0.78466 



Table A. 2. Correlation coefficients Kij in the error for- 
mula Eq. (A.2) 
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Kij 
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Kij 
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0.00273 
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1.00646 
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0.00000 
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-0.61883 
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0.00000 
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2.99961 
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4 


0.00000 
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-5.27980 
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1.66580 
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4 


12.58549 
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